PCT/EP20047003320 / 2003«>0596^wp^ 

Description V/ y ^ Rec'dPCT/PTO 2 4 OCT 200J' 

Method, configuration and c^puter program product comprising 
program code means and computer program product for 
determination of a selected position of a mobile 
communications device in a communications network 

The invention relates to determining a position of a mobile 
communications device in a communications network 
(localization) . 

With the increasing spread of mobile communications, the 
demand for additional services in the mobile radio system is 
also growing. 

Location Based Services in this case are taken to mean 
additional services of mobile radio operators which are 
offered or provided to users of the mobile radio system on a 
location basis, i.e. depending on a position or a geographical 
location of relevant users, for example location or distance- 
dependent usage tariffs or orientation aids for search and 
rescue services. 

Fundamental to a "Location Based Service" is therefore the 
localization or position determination of the user concerned 
or of their mobile communications device. 

Different techniques are known for such localization of mobile 
communications devices in communications networks, for example 
determining the position on the basis of delay time 
determination or delay time measurements of communications 
signals of a mobile wireless communications device to a base 
station of a communications network ([1], [2]) or a 
localization by means of a satellite-assisted system such as 
GPS. 
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The delay time-based position determination known from [2] is 
executed for a mobile telephone, generally a mobile station, 
in a GSM communications network (= Global System for Mobile 
Communications) ([7], [8], [9]) in accordance with a TDMA 
(Time Division Multiple Access) mobile radio technology. 

An individual mobile station, which has registered with a 
fixed local base station (a base station which is controlling 
the call) is currently assigned a free time slot in a TDMA 
frame by the base station. 

The communications signals intended for the mobile station 
concerned go to this time slot in signal packets, known as 
bursts, with a length of 15/26ms, from the base station or the 
communications signals or bursts sent from the mobile station 
must arrive at the base station. 

As a result of scattering the communications signals emitted 
by the base station find their way to the mobile station via 
different paths (multipath propagation) , in which case they 
are attenuated depending on their frequency. 

A receive field strength of the communications signals 
received by the mobile station is thus not only dependent on 
the distance from the mobile station to the base station, but 
also on the frequency and topographical circumstances between 
mobile station and base station. Therefore the individual data 
packets are sent on different carrier frequencies, which 
allows selective interference on one frequency to be 
distributed between a number of users. 

However a precise synchronization between mobile station and 
base station is needed to enable this to be done. This 
synchronization is additionally made more difficult by a 
mobility of the user, since the mobile station is now located 
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at changing distances from the base station and its 
communications signals exhibit different delay times. 

To equalize the different delay times and enable frame- 
synchronized data to be delivered to the base station, the 
mobile station measures the signal delay time to the base 
station and uses this to correct the start of sending its 
bursts . 

The signal delay time is encoded in a "timing advance" (TA) 
and exhibits a dependency on the distance between mobile 
station and base station controlling the call. 

64 stages are available for the TA which are (bit) encoded 
with the values 0 to 63 and represent the delay time. 

Since positions of base stations are known the position of the 
mobile station can be concluded from a TA or from the signal 
delay time. 

A measurement accuracy for determining the delay time relates 
to a bit duration, which is 48/13 in GSM corresponding to a 
single path length of approx. 554m. 

Determining the position of a mobile communications device in 
a UMTS (= Universal Mobile Telecommunication System) network 
is known from [3] . 

In the corresponding UMTS mobile radio standard on which the 
UMTS network is based, determining the position of a mobile 
radio device is already explicitly included in the standard or 
is required by it (TS 25.305 V3 . 1 . 0 : stage 2 "Functional 
Specification of Location Services in UTRAN" (release 99) , 
3GPP TSG-RAN-WG2, 2000) . 



Further methods for localization of a mobile communications 
device in a communications network are known from [4], [5] and 
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[6] . 



A non-linear filter, a Prior Density Splitting Mixture 
Estimator (PDSME) is known from [10] . 

This non-linear filter, this PDSME, is based on new Gaussian 
mix filter algorithms for non-linear system state equations 
and/or non-linear measurement equations. 

Usually linearization errors occur if non-linear system state 
equations and/or non-linear measurement equations are replaced 
by linearizations during the use of linear filter techniques. 

The PDSME reduces these linearization errors by splitting the 
a-priori densities into a number of Gaussian mix components 
each with smaller covariances. 

The PDSME can be applied both what is known as a "Prediction 
Step" for the non-linear systems state equations and also for 
what is known as a "filter step" for the non-linear 
measurement equations . 

Furthermore a measurement is known for the linearization error 
from [10] which is comparable with a Kullback-Leibler 
distance. In addition a desired estimation quality can be set 
for the PDSME. By limiting the maximum number all mix 
components for the PDMSE an upper limit can be specified for 
the computing effort. 

The localization methods mentioned have the disadvantages that 
the positions of the mobile communications devices determined 
by them are inexact and consequently affected by large 
uncertainties. More precise methods however require expensive 
additional devices and costly modifications to the 
communications network or networks and communications devices. 



The object of the invention is thus to make possible a 
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localization of a mobile communications device in a 
communications network which is accurate and is affected by 
the fewest possible uncertainties, which can be implemented as 
simply as possible and at low cost. 

This object is achieved by the method and the configuration as 
well as by the computer program with program code means and 
the computer program product for determining a selected 
position of a mobile communications device in a communications 
network with the features in accordance with the relevant 
independent patent claim. 

With the method for determining a selected position of a 
mobile communications device in a communications network with 
at least one first base station set up for a first 
communication with the mobile communications device by means 
of first communications signals, a non-linear communications 
model with a deterministic component and stochastic component 
is used. 

The non-linear communications model includes a deterministic 

and a stochastic component, 

the deterministic component of which describes a dependency 
between communications signals of a base station and 
positions of the mobile communications device and 
the stochastic component of which describes the uncertainty 
of the deterministic component. 

By using the first communications signal belonging to the 
selected position of the first base station and the non — 
linear communications model, the selected position is 
determined, with a non-linear Bayesian filter technique being 
applied to determine the position. 



By applying the non-linear Bayesian filter technique to the 
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communications model or to the communications signal, a 
possible location area of the selected position is determined 
from which the selected position is finally determined. 

The non-linear Bayesian filter technique in the invention is 

generally taken to mean the technique described below: 

The Bayesian filter technique corresponds to a continuation 
or modification over time (also computing the change over 
time) of known probability distributions of the system 
states and system outputs by the presence of new measured 
values . 

The inventive localization method is based on the idea of 
obtaining from communications signals between at least one 
base station (transmitter) and one mobile station (receiver) 
available in communications networks such as for example WLAN, 
GSM, DECT networks, distance-relevant parameters and from 
these geographical information, in this case a possible 
location or distance area or area of uncertainty of the mobile 
station. 

It should be noted that technologically the transmitter and 
receiver of the communications signals can also be implemented 
the other way round. The inventive mode of operation is not 
affected by this. In this case the base station would be the 
receiver and the mobile station the transmitter of the 
communications signals discussed. 

The inventive localization method is a further based on the 
knowledge that the communications signals received by the 
mobile station are dependent on the geographical position of 
the mobile station in the communications network or in 
relation to the base station. 

This dependency between at the received communications signals 
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or between a signal strength of the communications signals 
received (by the base station) and the geographical position 
of the mobile station can be described by what are referred to 
as measurement equations, i.e. by the communications model. 

These measurement equations or the communications model 
include a deterministic component, said deterministic 
component describing a dependency between the received 
communications signal and the position of the mobile station, 
and a stochastic component, said stochastic component 
describing an uncertainty of the deterministic component. 

This uncertainty can for example relate to the communications 
signal and/or to the dependency mentioned (model uncertainty) . 

The non-linear Bayesian filter technique applied in accordance 
with the invention to these measurement equations or the 
communications model for determining the possible location 
area of the mobile station or the selected position 
significantly improves the quality of the localization 
compared to usual filter techniques such as an Extended 
Kalmann Filter (EKF) . 

The non-linear Bayesian filter technique can be undertaken 
once or iteratively a number of times in a "filter step". 

In this possible location area or area of uncertainty, 
described using stochastic variables, such as a mean and/or a 
variance, the mobile station or its position is finally 
presumed. 

The presumption can be formulated mathematically by a 
characteristic value of the possible location area or area of 
uncertainty such as a focal point or an expected value. This 
can then be used as an estimation for the selected position of 
the mobile station. 
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A particular advantage of the invention lies in the fact that 
the localization is conducted on the basis of communications 
signals and known positions of base stations which occur in a 
normal operation of a mobile radio system and are available 
there. This enables expensive modifications and expansions as 
well as additional measurements of existing mobile radio 
systems or for existing mobile radio systems to be dispensed 
with. 

The computer program with program code means is set-up to 
execute all steps in accordance with the inventive method for 
determining a position, i.e. the inventive localization method 
if the program is executed on a computer. 

The computer program product with program code means stored on 
a machine-readable medium is set-up to execute all the steps 
in accordance with the inventive localization method of when 
the program is executed on a computer. 

The configuration and the computer program with program code 
means set up to execute all steps in accordance with the 
inventive localization method of when the programme is 
executed on a computer, as well as the computer program 
product with program code means stored on a machine-readable 
medium set up to execute all steps in accordance with the 
inventive localization method when the program is executed on 
a computer are especially suitable for executing the inventive 
localization method or one of its further developments 
explained below. 

Preferred developments of the invention are produced by the 
dependent claims. 

The developments described below relate to both the method and 
to the configuration. 
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The invention and the developments described below can be 
implemented both in software and also in hardware, for example 
by using a specific electrical circuit. 

Further the realization of the invention or of a development 
described below is possible through a computer-readable 
storage medium on which a computer program product with 
program code means is stored which executes the invention or 
development. 

Also the invention or any development of it described below 
can be realized by a computer program product which features a 
storage medium on which a computer program product with 
program code means is stored which executes the invention or 
development . 

For communication in a communications network between a mobile 
communications device (mobile station) , for example a mobile 
telephone and a base station, for example a circular antenna 
or a circular radiator or one or more sectoral antennas, data, 
the communications signals, is transmitted in signal packets 
known as bursts. 

Based on or using at the transmitted communications signals or 
signal packets various parameters relevant to distance can be 
determined which in their turn can be included as a basis for 
determining the possible location areas or distance areas . 

Such a distance-relevant, i.e. distance-dependent parameter is 
for example a field strength of a signal packet. 

The field strength exhibits a natural dependency on the 
distance between the mobile station and the base station 
(controlling a call) and consequently delivers information 
about the possible location area or distance area (area of 
uncertainty) of the mobile station. 
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This dependency between field strength and distance can be 
described by physical models which describe a propagation 
behavior of signals. 

For the determination of the selected position a 
characteristic value of the possible location area or area of 
uncertainty can be determined, such as a focal point or an 
expected value which is then used as an estimation for the 
selected position of the mobile station. 

The quality of the localization can be further improved by 
using a user model which describes a movement of the mobile 
communications device. Thus such a uSer model can limit the 
maximum step length or movement distance for a prespecified 
time step . 

The invention is especially suitable for use in at the 
environment of a digital cellular mobile radio system such as 
a GSM network, and there for example for localization of a GSM 

telephone (mobile telephone) . 

In this case when the invention is used only the data 
available to a mobile telephone will be used, in which case 
costly changes to not have to be made either to the GSM 
network or to the mobile stations in the GSM network. 

For example the positions of the individual base stations and 
their antennas as well as their characteristics which provide 
information about the coverage area of the relevant antenna 
are known from a GSM network. 

The mobile telephone for its part always stays in contact with 
the antennas that can be received for a correct connection 
set-up so that it can be assigned by the network the antenna 
which is most suited for a call. To this end it measures the 
field strengths of the receivable antennas as well as specific 
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signal delay times which are then also known. 

On the basis of this available information the mobile 
telephone is then localized in accordance with the inventive 
method of operation. 

The invention is also suitable for use in the environment of 
other digital cellular mobile radio systems such as a WLAN or 
of a DECT network (11], and there for example for localization 
of a DECT mobile telephone. 

The figures show an exemplary embodiment of the invention 
which is explained in more detail below. The exemplary 
embodiment is in this case subdivided into a basic part which 
presents the fundamentals of the inventive method of operation 
and a related part which clarifies or specifies on the basis 
of a concrete numerical example the inventive method of 
operation as well as the results. 

The Figures show 

Figure 1 a drawing in which it an assumed linear approximation 
for the non — linear communications model is 
presented; 

Figure 2 a drawing in which the uncertainty of the 

communications signal through noise is presented; 

Figure 3 a drawing in which an overview of the PDSME algorithm 
with a linearized filter step (upper section) and a 
prediction step (lower section) are shown; 

Figures 4a and b drawings in which an evaluation scenario 
for localization in accordance with the exemplary 
embodiment is shown; 



Figures 5a and b 



drawings in which a localization using an 



PCT/EP20047003320 / 2003P05965WOUS 



12 



extended Kalman filter is shown (first filter step 
Fig. 5a; 170th filter step Fig the. 5b); 



Figures 



6a and b drawings in which a localization using the 
PDSME in accordance with the exemplary embodiment is 
shown (first filter step Fig. 6a; 170th filter step, 
Fig. 6b; 



Figures 



7a and b drawings in which results of the EKF 

filtering and the PDSME filtering in accordance with 
the exemplary embodiment in relation to the correct 



results are shown; 



Figure 8 a drawing in which a number of Gaussian mix 

components which depend on filter steps are shown. 

Exemplary embodiment Localization of a DECT mobile telephone 
in a DECT network based on a non-linear filter technique 



The application of a non-linear Bayesian filter technique in 
the localization of a mobile radio telephone is described 
below. The application is described using as an example a 
localization of a DECT mobile telephone in a DECT network with 
a number of base stations. 

Signal strengths (as field strengths of the communications 
signals received by a mobile telephone and measured there (of 
a base station transmitting the signals) are dependent on the 
position of the mobile telephone in relation to the base 
station sending out the communications signals. This 
relationship is described through non-linear measurement 
equations (non-linear communications model) . 



Basics 



These non-linear measurement equations or this non-linear 
communications model comprises a deterministic component which 
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describes the signal strengths received as a function of the 
position and a stochastic component which takes account of 
model errors and measurement noise. 

In addition user models are also considered which bring with 
them knowledge about a (spatial) movement of users of the 
mobile telephones • 

A new non-linear filter technique, known as a "Prior Density 
Splitting Mixture Estimator" (PDSME) , assignable as a Gaussian 
mix filter algorithm significantly improves a localization 
quality compared to standard filter techniques such as an 
"Extended Kalman Filter" (EKF) (cf. Figures 5 to 7). 

Typical applications for localization tasks are determining 
the positions of mobile communications devices in WLANs, GSM 
networks or DECT networks. 

For position determination information should be able to be 
used in such cases which is already available during of a 
normal operation phase, i.e. for "normal" communication 
between mobile device and base station. 

The localization method known from [12] for mobile telephones 
in GSM networks is based on location-dependent signal 
strengths of the communications signals sent out by base 
stations, able to be received and measured by the mobile 
telephones. Because of the location dependency these are 
characteristic for the position of the signal recipient, in 
this case for the position of the mobile telephone or of a 
user of the mobile telephone 

During the normal operation phase the received signal 
strengths of all receivable communications signals of the 
relevant transmitter (base stations) are measured and compared 
by the mobile telephone. This is how a handover of the mobile 
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telephone between different transmitters or base stations is 
controlled. 

Based on these measurements made during normal operation 
localization methods can be developed in order to additionally 

determine the position of the mobile telephone in relation to 
the transmitter or in the GSM network. 

Further procedures are known from [13] and [17] for 
localization of mobile communications devices. These 
procedures for localization are based on an application of 
non-linear status estimations. 

The application of a new non-linear filter technique for the 
localization of DECT mobile telephones will be described 
below. 

To make it possible to localize mobile telephones, as is known 
from [18], a stochastic approximation model for the 
propagation of radio waves or receivable signal strengths of 
communications signals in a localization environment is 
identified or developed. 

A corresponding propagation model can also basically be 
obtained through physical propagation models of 
electromagnetic waves. 

An adaptation of purely physical propagation models to real 
(localization) environments is however very complicated since 
reflections or interference or other physical phenomena 
influence and falsify the assumed physical propagation. 

Since real propagation parameters which describes the 
electromagnetic (propagation) characteristics in the real 
localization environment in which the mobile telephone is to 
be localized are only partly known or determinable the use of 
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such physically-based propagation models for such real 
applications is excluded. 

Therefore a propagation model which is based on actual 
measurements of the (logarithmized) received signal strengths 
of the communications signals in the actual localization 
environment or which can be derived from it is to be preferred 
over theoretical physical propagation models. 

The development of such a measurement model or communications 
or propagation model derived from actual measurements requires 
a calibration within the framework of which parameters of the 
measurement model are defined. 

Calibration ("model generation phase") 

For such a calibration the receiver (used in the sense of a 
measurement device) i.e. the mobile telephone, is moved to a 
number of positions within the localization environment, for 
example at grid points obtainable by placing a grid over the 
localization environment. 

At each grid point the field strengths of all base stations 
which can be received at this point are measured and stored in 
a map of the localization environment- In this way the field 
strength distribution of each base station over the 
localization environment is determined. The field strength 
distributions can be modelled by the measurement model. 

Based on this model information (and on field strengths 
measured at an actual position) the actual position assumed by 
a receiver or by the mobile telephone in the localization 
environment can then be determined (localization in a 
localization phase) . 



Furthermore a description is given below of how a stochastic 
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measurement model can be identified in the model generation 
phase . 

This stochastic measurement model consists of an analytical, 
deterministic measurement function which describes the 
logarithmized received signal field strength as a function of 
position co-ordinates. In addition a stochastic component is 
identified which takes account of model uncertainty and also 
measurement noise . 

Such a stochastic approach to the localization takes account 
of the uncertainties in the position determined through 
probabilities or probability density functions. 

As well as the measurement model for the receivable 
communications signals the quality of the localization is 
further improved by the fact that a user model is included as 
well. Via this user model knowledge about a (spatial and/or 
temporal) movement of a user can be incorporated. 

Through such user models maximum "step lengths" of users in a 
time step which is defined via two associated or corresponding 
signal measurements (at the start and at the end of the times 
step) can be restricted. 

An additional improvement of the localization quality can be 
achieved by recursive position determination through a 
combination of multiple different consecutive measurements. 

The new stochastic approach used here for localization results 
in non-linear multi-dimensional measurement equations. An 
exact solution of the Bayesian filter problem needs to 
complicated non -Gaussian probability densities which describe 
the position sought along with its uncertainties. 

These probability density (functions) can be approximated by 



PCT/EP20047003320 / 2003P05965WOUS 

17 

Gaussian mix densities [19], [20], 

This approximation is undertaken here by a newly developed 
"Prior Density Splitting Mixture Estimator" (PDSME) which has 
been developed to make it possible for the user to set the 
quality of the localization - by contrast with usual, non- 
linear, known filter techniques such as an Extended Kalman 
filter (EKF) . 

However specifying an upper limit for a computation effort by 
restricting the number of Gaussian mix components cannot 
reduce the localization quality. In addition, by the 
representation of densities in the Gaussian mixtures a 
recursive method of operation both for the non-linear 
predictions steps and also for a non-linear filter steps is 
possible . 

Furthermore the first part describes how the stochastic 

localization problem can be formulated. The deterministic 
components and the stochastic components of the measurement 
model will be explained below. The use of Bayesian filter 
algorithms based on a Gaussian mix approximation for the exact 
probability densities used in the new PDSME will be described. 

In the second part results of a localization of a DECT mobile 
telephone in a DECT network will be presented and the quality 
of the localization using the PDMSE will be compared to that 
using a known EKF. 

Formulating the localization problem 

The localization task for a mobile communications device in a 
communications network can be divided up into two basic 
problems : 

1) the identification of the measurement model with a 
deterministic and a stochastic component in the model 
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generation phase; 
2) the determination of the probability density of the position 
sought in a Bayesian filter step in the localization phase. 

The measurement model 



Xk = 













Yk,N 






+ 





(1) 



describes the measurement of the logarithmized received signal 
strengths 



Yk, i; = 10 • logic 



ImW 



of the ^th transmitter, ^ = 1, ...^N, as a deterministic non- 
linear function h < (xk) of the position coordinates xk=[xl,k, 
x2,k] for each of the N different transmitters. 

The stochastic component vk of the measurement model is 
represented by an additive uncertainty. Additive uncertainties 
in the measurement model of the logarithmized receive signal 
strengths correspond to multiplicative uncertainties of the 
noise signal power P which are to be attributed to an 
influence of unspecified attenuations. 

Each transmitter can be identified by its own ID which is 
transmitted by the base station during a communication. 

This makes the localization problem considerably simpler since 
each measured value can be directly assigned to the 
appropriate non linear measurement equation. 

In the localization phase the non-linear filter algorithms 
will be applied in each time step k. However without the 
additive uncertainties vk the implicit, non-linear 
relationships between the vector with the measured 



PCT/EP20047003320 / 2003P05965WOUS 

19 

logarithmized received signal strengths yk of the transmitter 
can only be resolved numerically for example by a recursive 
least squares method. If the uncertainties are also taken into 
account the measurement equations can be used to update the 
position sought in the filter step. 

The further described method of operation for the localization 
is based on an approximation of the exact solution of the 
Bayesian filter step by Gaussian mixture densities 



with L components, defined by non-negative weights or expected 
values ]i and covariances C. 

If the measured value y is available the exact solution of the 
Bayesian filter step is produced: 



with the non-linear measurement function h{x), the additive 
uncertainty density functions f (v) and a normalization 
constant C. 

A- priori knowledge about the position sought is represented 



For each new measurement the previous a-posteriori density 
f(x|y) is interpreted as a new a-priori density so that the 
position sought or determined can be recursively updated. 

Similarly to the Bayesian filter step a prediction step which 
is used to describe the non-linear user model 




by the density function f(x). 



2lc+l ~ ^kO^k) i?k 



(2) 
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is likewise executed for approximation of the exact densities 
by Gaussian mixtures. 

The non^linear function a{x) is a deterministic model of a 
user movement. Uncertainties are again taken into account by 
additive noise. 

Stochastic modeling 

The measurement model is identified using measurements of the 
logarithmzed signal strengths for each transmitter at grid 
points of a grid which covers the localization environment. 
The deterministic (h(x))and the stochastic component (v) of 
the measurement model (1) are identified before the 
localization. 

The deterministic, analytical component h(x) is a measurement 
equation which describes the received signal power as a 
function of the position coordinates. The stochastic component 
V is a model for the uncertainty of the deterministic 
component. These uncertainties are made up of a spatial 
uncertainties corresponding to model errors and measurement 
noise over time. 

Deterministic measurement model of the logarithmized received 
signal strengths 

Fig. 1 shows that the acceptance of the logarithmized receive 
signal strengths over a distance of a number of meters is 
almost linear. In other documents similar assumptions in 
relation to the measurement model are referred to as "linear- 
loss-model" [24] or "linear-scope-model [18]]. 

Mathematically linear acceptance of the logarithmized receive 
signal strengths for two-dimensional position coordinates x 
can be described by N independent measurement equations. 
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The parameters mC and PC of the positive semi-definite 
quadratic form and the additive offset A< are identified for 
each transmitter ^=1, ...^N. 

The 6N parameters m^, P< and are determined by a method of 
recursive least squares in which the deviation between the 
approximated measurement model and the measured logarithmized 
signal strengths at the grid points measured during the model 
generation phase are minimized. To reduce the measurement 
noise over time the mean value over a number of measurements 
at each grid point Is used to calculate the average of the 
measured values. 

Stochastic modeling of the uncertainties 

The stochastic uncertainty model takes account of both the 
deviation between the approximated model described above and 
the true distribution of the logarithmized receive signal 
strengths measured at the grid points and the measurement 
noise over time. 

Model ucertainties 

In this localization approach the deviation between true 
logarithmic receive signal power and the deterministic 
measurement model h(Xk) is described for each transmitter by a 
Gaussian noise density with a mean value of //^J^ and standard 

deviation a^J^ . This Gaussian uncertainty is a representation of 

the approximation error of the deterministic components of the 
measurement model over the entire localization environment. It 
first represents the mean deviation over the localization 
environment between the model and the true receive power on 
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the basis of the incorrect assumptions for the deterministic 
part of the measurement equation. The uncertainties produced 
from large deviations between the true receive signal power 
and the measurement model are shown in Fig. 1. Secondly 
measurements have also shown that local deviations exist which 
could have been caused by reflections, non-homogeneous 
propagation of the radio waves and interference. These spatial 
variations of the logarithmic receive signal power can be seen 
from Fig. 2 in which measurements have been collected at a 
spacing of 2 cm. 

Measurement noise 

In addition to spatial variations of the receipt signal power 
measurement noise over time can be determined by evaluating a 
series of several different measurements at a fixed point. 
This measurement noise over time is also shown in Fig. 2 for 
three different measurements of the logarithmic receive signal 
power at each measurement point. This noise is in its turn 
approximated by a Gaussian density with a mean of fi^^^ and 

standard deviation al^^ . 

Combination of modeling certainties and measurement noise 

To obtain a simple model with the two uncertainties described 
above it is further assumed that the two uncertainties are 
independent. They can thus be modelled by a single Gaussian 
density fy^g(u^) for each transmitter 

^= 1,...,N, which is defined by the mean 
and the standard deviation 
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This model implies not only independent uncertainties but it 
is also assumed that the two uncertainties can be described 
without taking into account any positional dependency. 

Filter algorithm 

This section gives a brief overview of the PDSME (Prior 
Density Splitting Mixture Estimator) used for the localization 
of radio communications devices. In addition an adaptation of 
the measurement equations derived in Section 3 h(Xk) is 
introduced to simplify the computation or the PDSME estimator. 
Furthermore a prediction step for a simple user model is 
presented. 

The PDSME measurement updating 

The PDSME algorithm for the localization of radio 
communications devices shown in this step is based on the 
calculation of a ^linearized measurement update for Gaussian 
mixture densities. The measurement updating step of this 
filter algorithm is shown in a block diagram in the upper part 
of Fig. 3. Splitting is based on the calculation of the 
linearization error 



for each component of the Gaussian a-posteriori mixture 
density. This criterion is very similar to the Kullback- 
Leibler distance [25] 




(4) 
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and its 



at 



the mean of the ith component of the a-priori-density 
function. By calculating the linearization error (4) the 
Gaussian a-priori mixture components which contribute most to 
the approximation error of the a-posteriori density are 
identified. To reduce this linearization error these Gaussian 
mixture components are replaced by being split up into a 
number of mixture components with smaller covariances using 
splitting libraries which have been optimized offline (see 



After this "analytical resampling" of the a-priori densities, 
the filter step can be calculated by a bank of EKF, 
linearization of the measurement equation at the mean of each 
component of the representation of the Gaussian mixture of the 
a-priori density. 

Thereafter redundancy is reduced by a merging step by 
combining a number of Gaussian mixture components into a 
single Gaussian with negligible approximation error. 

Adaptation of the measurement model for simplifying the 
linearization error criterion 

For Gaussian measurement noise the linearization error (4) can 
be calculated analytically as a linear combination of moments 

of the densities fefcljffii|<) polynomial measurement 

equations h^ (Xk) . In this subsection an adaptation of 
measurement equation (3) is thus derived. After a few 
algebraic transformations the quadratic measurement equation 



Fig. 4) . 
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can be rewritten with a polynomial function with a modified 
"measured value" Zk/ < and the transformed uncertainty U^. 
Because of the non-linear transformation of the random 



:io„ /".f 



variable u^ the probability density function ^ ^ h/ no 
longer Gaussian. In the localization experiment in Section 5 
the exact moments of first and second order 0^ will be 

fuCivt) 

calculated to determine a Gaussian approximation of " ^ s/ , 
It should be noted that the moments depend on the measured 
value yk/c- They must thus be recalculated for each new 
measurement of the logarithmic receive signal power and are 
not unchangeable over time like the parameters described in 
Section 3 liu/c and Ou/^. 

User modeling by PDSME prediction step 

Like the filter step the PDSME can also be applied to the 
calculation of non-linear prediction steps. In a similar way 
to non-linear filter steps the calculation of an approximated 
prediction step also consists of an evaluation of a 
linearization error, a bank of linearized prediction steps and 
the reduction of the number of Gaussian mixture components in 
a merging step (see lower part of Fig. 3). In this work only a 
linear user model is taken into consideration. The predictions 
step can thus be calculated analytically since a-posteriori 
density has been approximated by Gaussian mixture density in 
the filter step. 

The prediction model consists of the linear status equation 
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with the mean ji^^ of the additive system noise co^ representing 
knowledge about possible directions and average step lengths 
of the movement of the user. The covariance matrix Cq of cok 
specifies an estimation for the distribution of the step 
lengths of the user. For each component i = 1,...,L of the 
Gaussian mixture density the predicted Gaussian mixture 
component is then described by the mean 



The estimated position can then be calculated as weighted 
overlaying 



of the mean values of all Gaussian mixture components. 
Part 2: Localization experiment 

In this section a real localization experiment for the 
validation of the localization approach described is presented 
for DECT mobile telephones. In these subsections the superior 
performance of the PDSME by comparison with the EKF (Extended 
Kalman Filter) , a widely-used standard approach for estimating 
the status of non-linear systems, is shown. 

Evaluation scenario 



and the covariance 




In this localization experiment N = 10 transmitters have been 
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placed in an internal space of approximately 30 m x 30 m on a 
floor of a building. In the model generation phase the 
logarithmic receive signal part of each transmitter has been 
measured on a grid with 1 m spacing between grid points. IN = 
10 measurement equations h4(Xk), consisting in total of 60 
parameters for the deterministic components and 20 parameters 
for the stochastic components jik have been identified. Apart 
from just a few areas which were influenced by the high 
attenuation of the radio waves as a result of steel and 
concrete walls, the assumed model is an appropriate 
approximation of the distribution of the receive signal power. 

The initial probability density of the position is selected as 
a Gaussian density with the initial mean value 

and initial covariance 

i.e. there is almost no prior knowledge available about the 
position. 

In the localization phase a measurement of each receivable 
transmitter is collected along a line, as shown in Fig. 4, 
with a spacing of Ar = 0,5 m between the true measurement 
points Pi, P2, Pi7- The position co-ordinates of the 
measurements have been determined to compare the basic 
actuality with the estimated positions. For the localization 
of the receiver the exact measurement positions have not been 
used either in the localization approach with the PDSME or 
with the EKF. After the measurement update for each receivable 
transmitter at a fixed point Pi,i = 1,...,17 has been 
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calculated, a prediction step in accordance with Subsection 
4.3 has been calculated. In the present example the simple 
user model is defined by the mean 

jJu = Om 

and the covariance 

i.e. there is no knowledge available about a preferred 
direction of the movement of the user. The covariance matrix 
only represents the distribution of the step lengths of the 
user. 

EKF filter (Extended Kalman Filter) 

When the EKF filter is used the measurement equation (5) is 
linearized at the mean of the a-priori density function. Fig, 
5 shows the results of the EKF for the first and the 170th 
filter step. For each filter step the measurement of the 
logarithmic receive signal power of a single transmitter is 
used for updating the estimated position. 

It can be seen that after the first filter step there is a 
significant estimation error. The true position of the 
receiver does not lie within the support values of the 
estimated a-posteriori density function, which is shown by the 
profile representation. Furthermore there is no crossing point 
between the true measurement equation and the estimated 
density. It is thus not possible to apply data validation 
processes to establish whether a measured value can be 
"explained" by the estimated density. 



PDSME (Prior Density Splitting Mixture Estimator) 



PCT/EP20047003320 / 2003P05965WOUS 

29 

Fig. 6 shows the a-posteriori densities for the first and 
170th filter steps estimated by the PDSME algorithm. The 
approximation of the non-Gaussian a-posteriori density in the 
first filter step which lies very close to the numerically 
calculated optimum Bayesian solution of the filter step is 
clearly much better than in the case of the EKF filter. Data 
validation procedures can thus be successfully applied. 

In Fig. 7 the expected values of the a-posteriori density 
functions calculated with the EKF and the PDSME are compared 
to the true position. If this figure is compared with Fig. 8 
it can be seen that for almost linear filter problems, i.e. if 
the covariance of the estimated position by comparison with 
the non-linearity of the measurement equation after a number 
of filter steps is small, EKF and PDSME yield almost the same 
results. In these cases the PDSME uses only a modest number of 
Gaussian mixture components whereas for strong non-linearities 
at the beginning of the localization experiment a higher 
number of approximation components is necessary to reduce the 
estimation error. The superior power of the PDSME filter by 
comparison with the EKF is also shown by the average 
estimation error 



over the Np = 170 filter steps which amount to 3,30 m for the 
EKF and 1,22 m for the PDSME. 

Conclusions 

In this work a stochastic approach to a localization of radio 

communications devices has been offered to which is based on 
the measurement of the logarithmic signal power of the 




/ 
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receivable transmitters by a mobile station. For each 
transmitter a stochastic measurement model has been identified 
which consists of a deterministic and a stochastic component. 
This measurement model has been used for the estimation of the 
position of the receiver through an innovative Gaussian 
mixture estimator which is based on the splitting of the a- 
priori density corresponding to a linearization error 
criterion. This criterion is very similar to the Kullback- 
Leibler distance between the true and the approximated a- 
posteriori density which is calculated by a linearization of 
the measurement equation. By applying this estimation 
procedure to the localization of DECT mobile telephones 
significant improvements to the estimation quality can be 
achieved if the PDSME is used instead of the standard 
approaches such as the EKF. Further improvements to the 
localization quality by identifying better deterministic 
measurement models and more exact identification of the 
measurement noise can be achieved which do not assume any 
independence between the different uncertainties mentioned in 
this work. 
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